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SUMMARY 

The  growth  and  collapse  of  the  turbulent  wake  behind  a  submarine  in  a 
stratified  ocean  is  studied  analytically  using  simple  phenomenological  models. 
The  present  analysis  closely  follows  the  previous  integral  approach  of  Ko, 
based  on  the  concept  of  local  similarity  and  turbulent  entrainment.  A  non- 
uniform  entrainment  rate,  obtained  by  considering  the  three  turbulent  velocity 
components  separately,  marks  the  main  difference  to  the  previous  approach. 
However,  the  present  analysis  is  still  limited  to  a  wake  without  swirl  and 
does  not  include  the  energy  leakage  through  the  internal  wave  field. 

A  series  of  parametric  studies  have  been  performed  to  study  the  effects 
of  various  operational  conditions  (such  as  stratification  and  speed)  as  well 
as  the  initial  conditions  on  the  wake  flow  field. 

The  solutions  were  found  to  be  sensitive  to  some  initial  conditions.  The 
physical  significance,  the  range  of  conceivable  values,  and  the  constraints  of 
these  initial  conditions  are  discussed.  The  differences  between  a  blunt  body 
wake  and  a  slender  body  wake  are  pointed  out,  based  on  these  parametric  studies. 

A  comparison  wi th  the  experiment  of  Schooley  and  Stewart  is  presented. 
Again,  an  extremely  high  initial  turbulent  intensity  and  degree  of  mixing  are 
needed  in  the  calculation.  In  addition  to  the  fact  that  the  body  is  practically 
a  blunt  one,  it  seems  to  suggest  that  an  excessive  amount  of  energy  (more  than 
required  to  balance  the  body  drag  for  maintaining  a  uniform  speed)  may  have  been 
imparted  to  the  wake.  A  much  different  result  is  expected  from  the  wake  of 
a  slender,  s ubma r i ne- 1 i ke  body. 
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1 •  Introduction 


Our  understanding  on  the  growth  and  collapse  of  the  turbulent  wake 
behind  a  submarine  moving  at  a  constant  speed  through  a  stratified  medium  has 
been  vastly  improved  in  recent  years.  The  simplified  integral  analysis  of  Ko 
(ref.  1)  has  provided  a  general  characterization  of  the  various  physical 
mechanisms  associated  with  this  complicated  flow  phenomenon.  However,  quite 
a  few  assumptions  and  approximations  in  that  simplified  analysis  were  under¬ 
stood  to  be  unsatisfactory  to  provide  a  reliable  prediction.  Therefore,  an 
improved  analysis  is  presented  here  with  the  hope  of  arriving  at  a  reliable 
first  order  prediction  tool  as  well  as  enhancing  our  physical  understanding  of 
the  submarine  wake. 

In  the  present  analysis,  the  background  ocean  is  assumed  to  be  "quiet" 
which  precludes  the  consideration  of  the  ocean  turbulent  diffusion  dominated 
far  wake  region,  (reference  2  gives  a  general  consideration  of  the  diffusion 
of  passive  scalars  behind  a  submarine  in  that  particular  region).  Furthermore, 
no  net  momentum,  either  in  the  form  of  a  swirl  or  a  lift,  is  assumed  to  be 
imparted  to  the  wake  In  the  present  analysis.  These  additional  submarine 
motion  associated  effects  will  be  considered  after  the  present  basic  model 
being  established. 
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2.  Technical  Approach  and  Formulation 


The  present  analysis  closely  follows  the  previous  integral  approach  of 
Ko  (ref.  1)  based  on  the  concept  of  local  similarity  and  turbulent  entrainment. 
The  main  advantage  of  the  integral  approach  over  a  finite  difference  approach 
is  the  simplicity  in  application  to  provide  a  quick  appraisal  of  the  various 
physical  mechanisms  associated  with  the  submarine  wakes.  The  shortcoming  of 
such  an  approach  is  the  inability  in  predicting  the  details  as  well  as  in 
giving  a  true  representation  of  the  region  very  close  to  the  body  where  a 
highly  non-similar  flow  field  is  expected. 

The  basic  flow  field  of  interest  is  shown  schematically  in  Figure  1.  The 
undisturbed  medium  outside  the  wake  is  assumed  to  be  horizontally  stratified 
with  the  density  given  by 


Po  =  P  (1-az) 


(1) 


with  a  being  a  constant, 
wake  is  taken  to  be 


The  corresponding  density  distribution  inside  the 


Pi  =  p  (1-gz) 


(2) 


with  £S  being  a  variable  equivalent  density  gradient.  In  the  integral  sense, 
the  wake  boundary  is  again  assumed  to  be  sharp  and  encloses  all  the  anamolies 
in  either  density,  turbulence  or  other  wake  identifying  quantities  (such  as 
dye).  The  possible  differences  between  a  dye  wake  and  others  will  not  be 
considered  here. 

It  is  again  assumed  that  there  is  no  direct  effect  of  turbulence  on  the 
inviscid  flow  field  in  the  cross  plane.  The  effect  of  turbulence  is  being 
accounted  for  through  the  moving  boundary  and  the  equivalent  local  density 
gradient.  Therefore,  as  derived  in  Reference  1,  with  the  locally  similar 
velocity  distribution 

v  =  yf^x) 


(3) 


w 


=  -zf^x) 


the  wake  boundary  can  be  shown  to  remain  elliptical  at  all  times  and  is  given 
by 


2  2 
2-  +  = 
2  v  2 

a  b 


1 


(4) 


- — - —  i  urt 


with  a  and  b  representing  the  major  and  the  minor  axis.  Furthermore, 
using  the  assumption  of  external  pressure  field  remaining  hydrostatic,  the 
cross  plane  momentum  equations  lead  to  a  governing  equation  for  the  "rate  of 
collapse"  f^ 


h2  df1  h2  h2  9 

Uo  (1  +  *2>  dx"  =  h  (a"6)  «  -  (1  -  \  )  f! 
a  a  a 


(5) 


The  concept  of  turbulent  entrainment  of  the  outer  fluid  leads  to  the 
following  governing  equations  for  the  rate  of  change  of  the  wake  boundaries 


U 


da 


,  =  af.  +  E  (x) 

o  dx  1  y 

db 


(6) 


b  ^  =  -bf.  +  E  (x) 
o  dx  1  z 


where  and  represent  the  entrainment  velocities  in  the  y  and  z 


directions  at  the  wake  boundary  y  =  a  and  z  =  b  .  The  previous  analysis 
assumes  a  more  or  less  uniform  entrainment  rate  by  linking  the  entrainment 
velocities  to  the  local  total  turbulent  intensity  level.  This  approximation 
is  definitely  undesirable  in  the  expectation  that  only  the  vertical  entrain¬ 
ment  rate  is  inhibited  by  the  effect  of  stratification.  Therefore,  an 
intrinsically  non-uniform  entrainment  rate  is  assume!  here  by  letting 


E  =  K.v 
y  1  m 


(7) 


E  =  K.w 
z  Ira 


with  and  w^  denoting  the  averaged  rms  turbulent  velocity  components  in  the 


y  and  z  directions.  Then,  the  rate  of  area  increase  is  given  by 


U  v'  w' 

-r-  4—  =  K.  <-S  +  -A 

A  dx  la  b 


(8) 
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Turbulent  Energy  Equations 

Instead  of  using  an  overall  turbulent  energy  equation,  three  separate 
energy  equations,  one  for  each  velocity  component,  are  used  to  account  for  the 
anisotropy  expected  in  a  stratified  flow.  The  governing  equations  with 
appropriate  simplifications  can  be  written  in  the  following  forms: 


ID  ,2  _ 
2  Dt  U 


— i — r  i — r  1  i  ' 

-uv  9y"uw  9z  -u9x 


±  D_  ,2 
2  Dt  V 


,2  3v  — — r  1  i 

-  v  —  -  v  w'  t - v  77—  - 

3y  3z  -  3y 


ID  ,2  — f— r  3w  ,  2  3w  1  ,  Sp*  '  g  -7— i 

2  dEw  ="vw  37  "w  77  -r"'  tr  -  E  -p  w 

p  p 

where  e'  denotes  the  viscous  dissipation  terms  and 


D  -3.-3j.~3 

—  =  U  —  +  V  T—  +  W  - — 

Dt  3x  3y  3z 


A  few  interesting  features  can  be  immediately  observed  from  equations  (9). 

(a)  Except  for  flow  with  an  appreciable  mean  swirl,  the  production  terms 

2 

effectively  appears  in  the  equation  for  u'  only.  The  exchange  of  energy  with 
the  mean  flow  is  mainly  carried  out  through  the  axial  component. 

(b)  The  pressure  velocity  correlation  terms,  which  serve  the  function  of 
redistributing  the  available  energy  among  the  three  components  in  an  attempt 

to  achieve  isotropy,  cancel  each  other  in  the  governing  equation  for  the  total 

turbulent  energy  _  _  _ 

e'  =  j  (u'2  +  v'2  +  w'2) 

(c)  The  effect  of  stratification,  with  the  Boussinesq  approximation, 
appears  in  the  governing  equation  for  w ' ^  only. 

The  modelling  of  each  term  is  briefly  discussed  in  the  following: 
Production 

The  simple  concept  of  an  eddy  viscosity  is  felt  to  be  adequate  for  the 
present  purpose.  It  is  assumed  that 
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The  eddy  viscosities  are  then  related  to  the  corresponding  turbulent  intensity 
and  integral  length  scale  by 


e 


y 


k'  X  (v'2)*5 
e  y 


e 

z 


K 


e 


X 

z 


(n: 


These  eddy  viscosities  are  further  assumed  to  be  constant  across  the  wake  but 

changing  in  the  axial  direction.  The  anisotropic  nature  of  a  stratified  flow 

is  built  into  Equation  (!1).  The  eddy  viscosity  constant  K  will  be  further 

e 

discussed  in  a  later  section. 


Dissipation 

The  modelling  of  the  energy  dissipation  term  follows  the  same  argument 
of  Townsend  (ref.  3)  but  equally  distributed  among  the  three  components. 
Specifically, 


then, 


1 

2  dt 


s'  * 


1.1  u 


,3 


=  0.3 


,3/2 


(12] 


Here  X  denotes  the  integral  length  scale  which  will  be  taken  to  be  the 
geometrical  average  of  X  and  X  i.e., 

y  z 


x 


(X  X  )** 
y  z 


(i3: 


Pressure  Velocity  Correlation 

The  pressure  velocity  correlation  terms  are  modelled  after  Rotta  (ref.  4) 
which  gives  - - — - — — 


where  denotes  an  universal  constant  having  a  value  of  about  ..6  as 
suggested  by  the  experimental  data  (refer  to  Hanjalic  and  Launder  (ref.  5) 
but  noted  that  some  algebraic  manipulations  are  required  to  arrive  at  the 
above  expression) . 

Stratification 

The  effect  of  stratification  is  also  included  using  a  turbulent  eddy 
diffusivity  model,  which  gives 

3pi 

-p7  ^  £p  jT  a 

Furthermore,  assuming  the  turbulent  Schmidt  number  in  the  vertical  direction 

to  be  unity,  we  have  _  i 

e  =  e  'v  X  (w'2)'* 
p  z  z 

Using  thece  simple  modellings,  the  three  governing  equations  (9)  can  be 
readily  integrat  id  across  the  wake.  Let's  define  the  area  avt raged  turbulent 
energy  in  the  three  principal  directions  by 


U'2 

A 

-  j  f 

U'2 

dA 

m 

2  A 

,2 

V 

A 

=  j  I 

v'2 

dA 

m 

2  A 

w'2 

A 

=  j  f 

„-2 

dA 

m 

2  A 

and  the  total  area  averaged  turbulent  energy  by 


e  =  u'2  +  v'2  +  w'2 

m  mm 


The  integral  governing  equations  become 


h  K2  «  ■  V  Dd2  +  T  > +  cp  t*  a(t  - 2  u:2) 


-  0.3A 


<L  fv.2 


fv'  A)  =  -2f  v 


,2A  +  C  ^  A  (v  -  2  v'2 


(18c) 


U  i-  (w'2  A)  =  2f  w'2  A  +  C  a  /2e  -  ,2, 

o  Jx  m  '  "lwm  A  +  cp  —  A  -  2w'  ) 


-  0.3a  — jp  -  f  Aw;  egb 

where  A  -  the  integral  length  scale  on  the  wake  axis.  A  =  c  (ab)^ 
in  the  present  model  with  C  =  0.125.  * 

cp  :  the  pressure-velocity  correlation  coefficient  which  is  related 
to  C  and  is  estimated  to  be  C  =08 

*  p  •  • 

kp  ;  \  =  °'18 

m 

Kp  :  an  integrated  constant  in  the  production  term  related  to  the  eddy 
diffusivity  constant  Kp  .  and  with  some  assistance  of  the  experi- 

mental  data  ,  we  have  K  =0  12 

P 

The  quantity  U,  represents  the  maximum  meat  velocity  excess  on  the  wake 
axis.  In  the  previous  analysis,  a  correlation  obtained  oy  Naudascher  (ref.  6) 
war  used  to  relate  this  quantity  to  the  local  turbulent  Intensity.  However, 
this  correlation  of  Naudascher  seems  to  be  unique  to  that  particular  experiment 
and  is  believed  to  be  a  weak  link  in  the  previous  analysis.  In  order  to  have  a 
better  assessment  of  the  energy  production  term,  it  is  decided  to  leave  U  as  a 
variable^  to  be  governed  by  the  mean  centerline  axial  momentum  equation.  d 
Let's  assume  the  mean  axial  velocity  field  to  be  given  by 

U  =  U  +  U ,  F  (£  — ) 

o  d  ’  b;  (19) 

The  mean  centerline  axial  momentum  equation  gives 

U  ~ -  =  -  (A  i2)  +  v  tt  /-Vm  ,  Wm,  „ 

0  dx  S>  3x  o  +  e  Ud  (X  +  T5  F  (20) 

where  F"  denotes  the  curvature  of  the  mean  velocity  profile  on  the  wake  axis 
and  is  the  eddy  viscosity  constant  defined  by 

Ey-Kev;a  ■  (2i, 
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Tlie  pressure  gradient  term  is  included  to  account  for  the  small  pressure 
gradient  field  induced  by  the  collapsing  motion,  which  c* n  be  easily  evaluated. 

For  the  case  of  a  momentumless  wake  and  zero  axial  pressure  gradient,  it 
can  be  easily  shown  that 

„  *2  ~ 

U^A  _  constant  (22) 


Using  eqn.  (8),  the  entrainment  constant  is  readily  given  by 

K.  =  0.5  K  F" 

1  e 


(23) 


Using  =  Kp  =  0.18  and  the  experimentally  measured  value  of  F"  a  25,  the 
entrainment  constant  can  be  estimated  to  be  Kj  =  2.25.  This  value  agrees 
fairly  well  with  the  growth  rate  obtained  in  Naudascher's  experiment. 


Mixing  Model 

Because  of  the  use  of  a  global  entrainment  concept  and  an  equivalent 
internal  density  gradient,  any  attempt  to  derive  a  mixing  equation  is  not  likely 
to  arrive  at  a  consistent  result  which  can  properly  account  for  the  effect  of 
turbulence  on  the  mixing  process  inside  the  wake.  For  that  reason,  we  have 
decided  to  use  a  phenomenological  model  to  obtain  a  governing  equation  for  the 
density  gradient  6  inside  the  wake. 

First  of  all,  the  total  potential  energy  inside  the  wake  is  given  by 


V  = 


upgU. 

___ 


Sa=s>-  ab3 

a 


(24) 


Therefore,  this  equivalent  linear  density  gradient  8  may  be  considered  as  an 
averaged  quantity  which  provides  a  measure  of  the  available  potential  energy 
inside  the  wake  at  any  given  station.  By  assuming  that  the  potential  energy 
per  unit  area  in  the  wake  is  directly  proportional  to  the  available  turbulent 
energy  per  unit  area  in  the  direction  against  the  gravitational  field,  the 
governing  equation  for  8  can  be  written  as 


w’ 

8b  =  ab  -  (a-B  )b  — y 
o  o  w 

mo 


(25) 
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It  should  be  noted  that  this  model  actually  is  not  applicable  to  any  ptemixed 
region,  such  as  the  laminar  experiment  of  Wu  (ref.  7). 

Nondimensionalization 

Following  the  same  nondimensionalization  procedure  of  reference  1,  the 
following  characteristic  quantities  are  used: 

Characteristic  length  L  =  initial  wale  size 

Characteristic  velocities  Uq  =  sub  speed 

^ref  =  eo  =  turbulent  intensity 

-1  -k 

Characteristic  time  N  =  (ag)  2 

Furthermore,  using  th>.  aster ik  to  denote  the  nondimensional  quantities,  the 
following  eight  dependent  variables  are  defined  by 

A  A 


a  = 

>  a/L  ,  b  =  b/L 

Wake  ciraensions 

a 

fl 

=  fj/N 

Collapsing  rate 

a 

u 

m 

-  u'/U  - 
m  ref 

A 

V 

ffi 

-  v' /U  f 
m  ref 

Turbulent  velocities 

A 

w 

m 

=  w'/U  f 
m  ref 

A 

e  = 

:  8L 

Internal  density  gradient 

U*  =  U^/U  Mean  centerline  axial  velocity  excess 

The  problem  is  then  reduced  to  finding  the  variation  of  these  unknowns,  which 
describe  the  wake  evolution,  as  functions  of  a  nondimensional  distance  defined 


where  =  »  the  Froude  number  based  on  the  body  diameter,  D  and  the  free- 

stream  velocity,  Uq.  These  unknown  quantities  are  governed  by  the  following 
set  of  first  order  ordinary  differential  equations: 

Entrainment 


da  *  *  * 

— r  =  a  f.  +  K.F  v 
,  *  1  1  t  m 

dx 


(27a) 


db  ,  *  *  * 

— t  =  -b  f  +  K.F  w 

dx  1  1  m 


( !/b) 


Cross  plane  momentum 


*2  df  *  *2  *2 

(1  +  J±.  .  a  .  JL  ,  ki  .  a  .  LI,  f  *2 

*2  dx*  «*  *2  *2  1 

a  a  a 


Turbulent  Energ 


*  *  a/2 

a*2**  ^  *  *2  v  w  *  *.  * 

d?  (um  a  b  }  =  Kp  “2  3  b  Ud  (-  +  7*>  -  °*3Ft  - pr- 
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(27c) 
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(27d) 


d  *  *  * 

*»  >  ' 
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(27f  ) 


Mixing 


d  (  *u\  a*  db 

“ *  k  b  >  “  “  — 

dx  dx 
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Mean  centerline  axial  momentum 


dU. 
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dx 
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It  is  clear  from  these  equations  that  two  nondimens ional  parameters  govern  the 
solutions:  the  turbulent  Froude  Number  F  defined  by 


F 

t 


ref 

NL 


(28a) 


and  the  initial  turbulent  intensity  parameter  given  by 


ref 

U 


(28b) 


The  turbulent  Froude  number  is  related  to  the  ordinary  Froude  number  by 


F 

t 


UN  (£>  Fr 


(28c) 


The  solutions  to  this  set  of  ordinary  different  equations  can  be  obta’ned 

*  * 

numerically  with  a  prescribed  set  of  initial  conditions  at  x  =  xq  .  In  fact, 
instead  of  appearing  explicity  in  the  governing  equations ,  a  few  fundamental 
parameters  to  this  problem  are  hidden  in  the  form  of  initial  conditions. 
Therefore,  a  full  discussion  of  the  physical  significance,  the  proper  choice  of 
values  and  the  oanstraints  on  the  initial  conditions  is  given  in  the  following: 

Initial  Conditions 
*  * 

(a)  a  ,  b  Initial  wake  dimension.  With  the  initial  wake  dimension  L 
o  o  -  *  * 

being  chosen  for  normalization,  both  a  and  b  can  be  set 
equal  to  one  if  the  wake  is  assumed  to  be  circular  initially. 
However,  experimental  observations  have  indicated  a  relatively 
strong  stratification  effect  on  the  initial  wake  aspect  ratio 

"k  k 

(b  /a  ),  the  effect  on  the  overall  solution  will  have  to  be 
o  o 

studied  parametrically. 


(27h) 


■«,  i * 


•4  V 
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('b)  f10  Initial  rate  of  collapse.  Assuming  that  the  fluid  particles 

do  not  possess  an  initial  velocity  field  in  che  cross  plane, 
the  initial  value  of  f^  can  be  readily  set  as  zero.  A  para 
metric  study  is  again  required  to  assess  its  significance  on 
the  solution. 


(c) 


mo 


mo 


w 


mo 


Initial  anisotropy.  Using  e  **  as  the  reference 


quantity,  it  requires  u  *^  +  v  +  w 

mo  mo  mo  o 

Thus,  the  relative  magnitude  of  these  quantities  represents  the 
initial  degree  of  anisotropy.  For  an  initially  isotropic 
turbulent  field,  *2  *2  i 


*2 


* 

e  =1 


mo 


=  v 


mo 


=  w 


mo 


(d)  Udo  Initial  mean  axial,  velocity  excess.  For  a  wake  of  zero  net 

axial  momentum,  this  value  represents  the  amount  of  mean 
energy  flux  remaining  in  the  wake  at  the  initial  station. 

(e)  Initial  mixing  effectiveness.  This  quantity  measures  the 

initially  available  potential  energy  in  the  wake.  A  combina- 
*  * 

tion  of  0o  and  f1Q  ,  gives  the  total  initial  energy  being 
imparted  to  the  wake  which  is  expected  to  control,  to  a  large 
extent,  the  collapse  process.  S.nce  little  is  known  about  the 
degree  of  initial  mixing,  it  is  decided  to  set  fj^*=0  in  the 
present  report  in  order  to  simplify  the  discussion. 

At  this  point,  we  seem  to  be  facing  a  large  number  of  possible  combinations 
of  parameters  and  the  proper  choice  of  values  becomes  extremely  critical  in 
assessing  the  potential  utility  of  such  an  analysis.  It  is  worthy  to  note  that 
thih  iS  exactly  ,;he  Primary  goal  of  this  type  of  simple  integral  approach  to 
bring  out  the  important  physical  variables  and  governing  parameters  in  a  complex 
flow.  A  sensitivity  study  of  these  parameters  must  be  performed  in  order  to 
identify  their  relative  significance  to  the  solutions.  Existing  experimental 
data  provide  a  reasonable  estimate  on  a  few  parameters,  but  additional  experi¬ 
mental  information  is  required  in  order  to  obtain  a  reliable  selection  of  the 
initial  conditions  for  future  application.  On  the  other  hand,  these  initial 
conditions  can  not  be  arbitrarily  assigned  even  though  the  governing  equations 
may  not  have  posed  any  restriction  on  the  values.  The  physical  constraints  on 
the  initial  conditions,  based  mainly  on  energy  consideration,  are  briefly 
discussed  in  the  following. 
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Constraints  on  the  initial  conditions 


(a)  Experimental  data  shows  that  at  an  initial  station  where  appreciable 


mean  shear  still  exists, 


max  ^  For  a  Gaussian  distribution,  the 


maximum  axial  turbulent  velocitv  intensity  can  be  related  to  the  area  averaged 


total  turbulent  intensity  by  e^  =  0.625  u'  Hence,  at  an  initial  station 

max 


not  too  far  downstream  from  the  body  where  the  momentumless  nature  becomes 
important,  the  initial  mean  axial  velocity  excess  is  related  to  the  initial 
turbulent  intensity  by 


U 


do 


-v  v  u 


N 


(29) 


with  K  being  about  2.  As  discussed  previously,  this  constant  K  may  be 
if  excessive  energy  still  remains  in  the  mean  flow  at  the  initial 
station.  To  some  extent,  this  "constant"  K  should  depend  on  the  configuration 
in  achieving  a  flow  with  no  net  axial  momentum.  For  example,  Naudascher's  type 
experiment  with  the  thrust  being  supplied  by  a  jet  is  expected  to  have  more 
residual  mean  energy  ns  compared  with  a  propeller  driven  body.  Additional  con¬ 
sideration  and  more  experimental  information  are  required  to  provide  a  better 
estimate  of  the  relative  energy  content  in  the  mean  and  in  the  turbulent  field. 
For  the  present  purpose,  the  crude  estimate  of  K  =  2  is  used  and  a  sensitivity 
study  will  be  performed  and  discussed  in  the  next  section. 

(b)  The  initial  total  turbulent  energy  flux  is  given  by 


E  =  pU  U  ‘ 
to  o  ref 


On  the  other  hand,  for  a  body  moving  at  an  uniform  speed  U  ,  the  total 


energy  imparted  to  the  wake  can  be  written  as 


Edrag  =  U0(DraS> 


I  pUc3  S*b 
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where  is  the  drag  coefficient  of  the  body  based  on  the  frontal  area  A^. 

Then  the  percentage  of  energy  goes  into  turbulence  is  given  by 

ft°  ,  1-  U2  (i)2  (30) 

Edrag  CD  8  ° 

A  typical  number  of  10  v20%  would  be  reasonable  for  a  slender  body  wake  without 
extraneous  energy  sources.  Therefore,  the  values  of  can  not  be  chosen 
arbitrarily.  It  is  also  interesting  to  note  that  an  appreciable  difference  in 
the  magniture  of  b’N  is  expected  between  a  slender  body  wake  and  a  blunt  body 
wake  because  of  the  difference  in  C D*  In  addition,  the  ratio  of  the  total 
turbulent  energy  to  the  total  energy  input  can  also  depend  on  the  body  config¬ 
uration  which  makes  the  choice  of  quite  different  for  blunt  and  slender  body 
wakes.  Further  discussions  will  be  given  in  the  next  section. 

(c)  The  total  initial  potential  energy  flux  can  be  written  as 

v  .!f£^Ab2 

o  8  a  o  o 

Then, 

_ L-  (V  T  '2  (1  -  —  )2  (31) 

Edrag  C/D  r 

Again,  this  ratio  poses  a  lower  bound  for  the  conceivable  values  of  (^/a).  We 
don't  really  have  a  good  estimate  on  the  percentage  of  available  energy 
appearing  in  the  form  of  a  potential  energy.  However,  it  must  satisfy  the 

inequality. 


<  1  - 


where  E  indicates  the  energy  associated  with  the  mean  flow  at  the  initial 
mean  . .. 

station.  It  may  also  be  noted  that  CD  again  plays  an  important  role  as  well  as 

the  Froude  number . 


3.  Results  and  Discussions 


3. 1  Nonstratified  Wakes 

As  a  check  for  the  empirical  constants  used  in  the  present  analysis, 
a  calculation  has  been  made  for  a  momentumless  wake  in  a  homogeneous  medium. 

The  initial  conditions  at  x/D  =  5  were  taken  from  a  recent  wind  tunnel  experi¬ 
ment  of  Gran  (ref.  8).  Since  the  data  is  still  being  reduced  at  the  present 
time,  no  specific  comparison  is  shown  here.  Figure  2  shows  the  growth  of  the 
wake  radius  as  a  function  of  (x/D).  Approximately  r  ^xU'  which  agrees  well 

with  the  preliminary  data  of  Gran  as  well  as  some  recent  FRI  towing  tank  data 
by  Fao  et  al.  (ref.  9) 

Figure  3  shows  the  decay  of  the  mean  centerline  axial  velocity  excess 

which  follows  a  nearly  x  law.  This  is  certainly  consistent  with  the  wake 

growth  rate  for  a  momentumless  wake.  Preliminary  data  of  Gran  shows  a  similar 

“2 

decay  while  a  much  faster  decay  rate  (wx  )  was  indicated  by  Naudascher.  It  is 

generally  agreed  that  the  mean  flow  measurements  of  Naudascher  are  not 

consistent  with  the  other  turbulent  measurements.  Therefore,  we  feel  that  the 

consistency  with  the  momentumless  nature  of  the  flow,  which  is  exhibited  by  the 

present  solution,  is  more  important  and  the  disagreement  with  Naudascher 's  mean 

flow  data  will  be  ignored.  Figure  4  gives  the  decay  of  the  turbulent  intensity 
*  -1 

um  ,  which  shows  a  nearly  x  dependence.  This  decay  rate  is  in  reasonable 
agreement  with  Naudascher  but  seems  to  be  slightly  faster  than  the  measurements 
of  Gran.  However,  as  noted  in  this  figure,  the  decay  rate  is  slowing  down  at 
increasing  x/D,  and  the  no-swirl  assumption  of  the  present  analysis  should 
account  for  part  of  the  faster  decay  of  the  turbulent  intensity. 

In  summary,  the  calculations  for  the  nonstratified  wake  has  provided 
enough  confidence  on  the  numerical  constants  used  in  the  present  analysis  and 
the  remaining  section  will  be  dealing  with  the  stratified  wakes  only  using 
these  prefixed  constants. 

3.2  Typical  Results 

Figure  5  shows  a  typical  calculated  variations  of  the  wake  height 
*  *  * 

(b  )  and  width  (a  )  as  a  function  of  the  nondimensional  distance  x  for 
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U„  =  e^/U  =0.1 
N  oo 

8  /a  =  0.5 

o 
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The  results  indicate  that  the  wake,  initially  circular,  grows  nearly 

axisymmetrically  before  reaching  a  maximum  height  and  then  "collapses".  The 

general  characters  of  the  solution  remains  the  same  as  the  results  of  the 

previous  analysis  (ref.  1).  The  phenomenon  of  "over  collapse"  seems  to  remain 

which  is  partially  caused  by  not  including  the  internal  wave  energy  losses. 

T.t  is  more  interesting  to  note  the  decay  of  the  three  turbulent  components 

Wi  n  is  shown  on  Figure  6.  The  calculations  were  performed  by  assuming  an 

equipartician  of  the  turbulent  energy  among  the  three  modes.  As  seen  in 

Figure  6,  because  of  the  mean  shear  production  term,  energy  is  fed  into  u* 

m 

initially.  JThen  the  stratification  effect  appears  and  starts  to  take  energy 
away  from  w^  .  The  tendency  towaid  isotropy  term  is  not  powerful  enough  to 
compensate  for  the  energy  gain  by  u  and  the  energy  drain  by  w*  .  Therefore, 
an  anisotropic  turbulence  field  is  developed.  In  fact,  at  soma  distance 
beyond  the  point  of  collapse,  the  turbulence  becomes  essentially  two  dimensional 
with  nearly  equal  amounts  of  energy  in  the  u*  and  the  v'  components.  In  the 
present  entrainment  model,  the  vertical  entrainment  is  practically  shut-off 
while  the  horizontal  entrainment  continues.  This  fact  is  reflected  in  Figure  7 

?Hk 

where  the  wake  area  is  plotted  against  x  .  In  the  previous  model,  an  asymptotic 
area  is  reached  at  a  short  distance  beyond  the  point  of  collapse. 

3.3  Effect  of  initial  anisotropy 

The  initial  turbulence  field  is  expected  to  possess  a  certain 

degree  of  anisotropy  either  as  a  result  of  the  stratification  effect  on  the 

body  boundary  layer  or  caused  by  the  propeller.  The  effect  of  these  initial 

anisotropy  on  the  flow  has  been  studied  by  setting  u*2:  v*2:  w*2  equal  to  1:1:1 

(the  typical  case  studied  above),  1:2:2  and  2:2:1.  Except  for  a  small  difference 
*  *  *  * 

ln  um»  vm  and  wm  at  small  x  ,  the  three  sets  of  solution  are  almost  indistin¬ 
guishable  from  the  results  shown  on  Figures  5-7.  Based  on  these  findings,  it  is 
concluded  that  the  initial  anisotropy  tends  to  be  quickly  removed  by  the 
tendency  toward  isotropy  terms.  Thus,  with  the  exception  of  a  strongly 
stratified  environment,  the  effect  of  anisotropy  on  the  wake  evolution  is  small 
and  an  equipartician  of  the  initial  turbulent  energy  among  three  modes  is 
assumed  throughout  the  remaining  calculations. 
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3.4  Effect  of  initial  degree  of  mixing 
As  stated  previously,  the  initial  density  gradient  g^  represents  a 

measure  of  the  total  amount  of  potential  energy  deposited  in  the  wake.  There¬ 
fore,  its  effect  on  the  wake  collapse  is  expected  to  be  quite  pronounced  as 

*  * 

demonstrated  on  Figure  8.  Three  sets  of  a  and  b  are  shown  on  Figure  8,  for 

6Q/a  =  0  (perfectly  mixed),  0.5,  and  0.9  (poorly  mixed).  The  results,  though 

expected,  are  strikingly  different.  In  particular,  for  gQ/a  =  0.9,  the 

"collapse"  is  hardly  "v  jible"  with  only  the  maximum  wake  height  being  limited. 

*  *  * 

By  defining  th=  distance  to  collapse  as  x  where  b  =  b  ,  it  is  easily 
^  c  max 

observed  that  x^  increases  rapidly  with  increasing  8Q/ct  .  This  type  of 

finding  has  been  extensively  used  recently  in  a  qualitative  manner  to  explain 

the  differences  in  the  various  experiments.  It  is  interesting  to  note  that  the 

effect  of  initial  mixing  on  the  overall  wake  evolution  is  quite  small  as  all 

three  sets  of  solution  curves  give  nearly  the  same  growth  rate  of  the  wake  area. 

The  variatic u  of  g/a  is  shown  on  Figure  9,  for  four  different  initial 

density  gradients.  It  should  be  cautioned  that  this  initial  mixing  parameter 

g  is  different  from  that  used  in  the  previous  analysis,  where  the  equivalent 
o 

density  gradient  was  kept  at  constant  before  collapse  as  shown  in  Figure  9  for 
the  purpose  of  comparison. 

The  intermediate  value  of  g  /a  =  0.5  has  been  used  for  most  of  the  para- 

o 

metric  studies  performed  here.  Some  definitive  measurements  for  the  mixing 
effectiveness  behind  the  full  scale  body  are  needed  to  provide  a  better  estimate 
of  the  values  to  be  used  for  prediction  purpose.  At  this  point,  we  can  only 
note  that  for  a  slender  body  with  =  0.2,  the  potential  energy  deposited  in 
the  wake  ranges  from  .05%  to  5%  of  ^(jrag  f°r  the  cases  shown  on  Figures  8  and  9. 

3.5  Effect  of  stratification  and  speed 

The  effects  of  stratification  and  speed  ran  be  grouped  into  a  single 
U 

parameter  Fr  =  »  the  Froude  number.  Figure  10  shows  the  wake  height 

variation  for  a  range  of  Froude  numbers.  For  fixed  values  of  and  gQ/a,  the 

maximum  wake  height  increases  with  increasing  F  (faster  body  speed  or  weaker 

r  * 

stratification).  The  nondimens ional  location  of  collapse,  x^  also  increases 

with  increasing  F  .  This  is  different  from  the  previous  analysis  which  gives 
* 

x  ;  constant  for  a  wide  range  of  Froude  number.  Furthermore,  the  degree  of 
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collapse  (measured  by  the  rate  of  collapse  at  the  maximum  height)  decreases  with 

£  Y  1 

increasing  F  .  Since  the  nondimensional  distance  x  =  — — ,  the  effect  is  much 
r  D  F 

more  evident  in  the  physical  coordinates  as  shown  on  Figure  11.  The  relative 

degree  of  collapse  can  be  readily  assessed  in  such  a  plot. 

The  explanation  for  the  strong  effect  of  Froude  number  can  be  obtained 

from  Figure  12  where  the  normalized  wake  area  is  plotted  against  x*  for  several 

Froude  numbers.  For  increasing  F  ,  the  relative  effect  of  stratification  is 

much  weaker,  the  vertical  entrainment  becomes  less  restricted  and  more  fluid  is 

being  entrained  into  the  wake.  The  effect  of  Froude  number  on  the  decay  of  the 

mean  centerline  axial  velocity  is  shown  on  Figure  13.  The  decay  of  U  tends 

d 

to  slow  down  for  increasing  stratification  (F  decreases).  This  general 
behavior  of  the  effect  of  stratification  on  the  mean  velocity  decay  has  been 
observed  experimentally  by  Pao  et  al.  (  ref.  9).  A  direct  comparison  will  be 
performed  in  the  future  when  more  experimental  data  becomes  available. 

3 . 6  Effect  of  initial  turbulent  intensity 

By  keeping  F^  equal  to  a  constant  (Fr  =  10),  Figure  14  shows  the 

effect  of  UN,  the  initial  turbulent  intensity  level,  on  the  wake  height.  The 

results  are  similar  to  the  previously  discussed  effect  of  Froude  number;  with 

increasing  UN,  Ffc  =  increases,  b^  increases,  x£  increases,  and  the  degree 

of  collapse  slightly  weakened.  However,  the  interpretation  is  entirely  different. 

First  of  all,  there  is  no  stretching  of  the  physical  distance  x  for  the  four 

diffe^6nt  values  of  U^.  Therefore,  their  effect  on  the  degree  of  collapse  is 

relatively  weaker.  Secondly,  as  discussed  in  the  previous  section,  the 

magnitude  of  is  a  measure  of  the  initial  total  turbulent  energy.  Also  shown 

on  the  curves  are  the  corresponding  values  of  /E,  for  each  U„.  The  bodv 

to  drag  N  J 

drag  coefficient  is  purposedly  left  unspecified.  For  a  slender  body,  the  drag 

coefficient  is  expected  to  be  about  0.1  m0.2.  If  the  ratio  E  /E,  is 

to  drag 

assumed  to  be  of  the  order  of  10  ^20%,  a  typical  value  of  UN  :0.05  would  have 
been  reasonable  for  a  slender  body  wake.  On  the  other  hand,  CD  ;  1  for  a  blunt 
body.  Thus,  a  higher  value  of  can  be  expected  and  its  effect  on  the  wake 
is  clearly  demonstrated  on  Figure  14.  However,  can  not  be  assigned 
arbitrarily  as  indicated  by  the  curve  for  UN  =  0.4.  Solution  can  still  b' 
obtained  but  more  energy  than  available  will  have  to  be  stored  in  the  turbulent 
field  even  for  a  blunt  body  wake.  This  is  certainly  unrealistic  unless  some 
extraneous  sources  of  energy  have  been  introduced  into  the  flow  field. 
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3.7  Effect  of  integral  scale  length 

Instead  of  using  a  separate  length  scale  equation  for  the  turbulent 
integral  scale,  an  algebraic  relationship  between  the  integral  length  scale 
and  the  equivalent  wake  radius  was  used,  i.e. 

A  =  C_^  (ab)^ 


After  searching  through  a  limited  amount  of  experimental  information,  the 

value  of  C.  was  taken  to  be  0.125.  The  effect  of  this  constant  on  the  wake 
x 

dimension  is  shown  on  Figure  15,  for  a  +  40%  variation  of  C^.  It  should  be 

noted  that  this  is  only  a  qualitative  assessment  of  the  effect  since  the 

integral  length  scale  only  affects  the  dissipation  and  the  pressure  velocity 

correlation  terms  in  this  model.  An  increasing  length  scale  implies  a  slower 

decay  as  shown  on  Figure  16  where  the  decay  of  turbulent  energy  level  is  plotted 
* 

verses  x  .  To  be  more  rigorous,  the  effect  of  integral  length  scale  on  the 
energy  production  term,  which  has  been  absorbed  in  the  constants  here,  must  be 
included.  This  modification  is  expected  to  accentuate  the  effect. 

It  is  also  interesting  to  note  that  the  experimental  observations  of  Gran 
seems  to  indicate  a  larger  for  a  drag  wake  than  a  momentumless  wake.  The 

corresponding  consequences  of  this  difference  may  be  inferred  from  these  results. 

3.8  Effect  of  initial  mean  velocity  excess 

* 

The  initial  value  of  gives  a  measure  of  the  amount  of  energy 

remained  in  the  mean  flow.  Since  the  transfer  of  energy  from  the  mean  to  the 

turbulent  field  is  relatively  inefficient  in  a  turbulent  wake  model,  the  effect 

on  the  wake  solution  is  quite  small  as  indicated  on  Figure  17  where  two  sets  of 

* 

curves  are  shown  bounding  four  calculations  for  =  0.05,  0.1,  0.2  and  0.3. 

The  corresponding  mean  velocity  decay  is  shown  on  Figure  18.  This  comparison 

* 

points  out  that  a  crude  estimate  on  the  initial  value  of  is  quite  adequate 

except  for  the  prediction  of  . 

3.9  Effect  of  initial  wake  geometry 

There  are  at  least  two  reasons  for  the  initial  wake  geometry  to  be 
non-circular  as  assumed  so  far.  First  of  all,  the  basic  body  geometry  is  not 
exactly  axisymmetric  and  small  angle  of  attack  can  have  a  strong  influence  on 
the  initial  wake  shape  for  a  slender  body.  Secondly,  the  effect  of  stratifica¬ 
tion  on  the  body  boundary  layer  is  expected  to  limit  the  vertical  growth.  Hence, 
a  larger  aspect  ratio,  a/b  is  expected  with  increasing  stratification.  This 
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dependence  of  the  initial  wake  geometry  on  stratification  has  been  observed 
qualitatively  by  Pao  et  al.  (ref.  9);  and  its  effect  on  subsequent  wake 
development  must  be  explored. 

Figure  19  shows  the  effect  of  the  initial  wake  geometry  on  the  downstream 
wake  evolutions.  The  initial  area  is  kept  at  constant  for  the  three  cases.  The 
effect  on  the  total  wake  area  variation  is  quite  small.  However,  the  effect  on 
the  apparant  wake  collapse  is  surprisingly  strong.  Practically  no  collapse 


in  the  classical  sense  can  be  stated  for  the  case  2.  On  the  other  hand,  a 
much  stronger  collapse  at  a  smaller  x  occurs  for  the  case  3  where  (b/a)  >1. 


This  result  provides  a  partial  explanation  to  the  much  weaker  collapse  being 
observed  recently  in  comparison  with  the  previous  experiment  of  Schooley  and 
Stewart  (ref.  10)  who  used  a  fairly  blunt  body,  for  which  a  more  or  less  circular 
initial  wake  shape  was  expected.  Additional  discussions  on  the  experiment  of 
Schooley  and  Stewart  will  be  given  in  the  following. 


3.10  Experiment  of  Schooley  and  Stewart 

In  the  previous  report  of  Ko,  a  comparison  with  the  laboratory 
experiment  of  Schooley  and  Stewart  was  presented.  Extraordinary  high  initial 


turbulent  intensity  (UN  =  0.45)  was  required  to  match  the  initial  wake  growth 


and  the  maximum  wake  height  observed  in  the  experiment.  It  is  worthy  to  note 
that  since  the  averaged  value  is  used  in  this  type  of  integral  analysis,  a  much 
higher  local  maximum  intensity  is  expected.  For  =  0.45,  with  a  Gaussian 
distribution,  we  have 


=  0.72 


a  nonrealistically  high  value.  Therefore,  an  attempt  to  use  the  present 
analysis,  together  with  the  understanding  achieved  through  the  paiametric 
studies,  to  compare  with  and  to  interpret  the  results  of  Schooley  and  Stewart 
is  given  here. 

The  flow  conditions  for  the  experiment  are: 


N  =  2.19  sec 


45  cm/sec 


L  =  3  cm 


F  =  -^  =  6.8 
r  NL 


A  larger  initial  wake  dimension  L  was  chosen  with  the  hope  of  getting  slightly 
farther  away  from  the  body  where  a  complicated  flow  field  is  expected  from  the 
relatively  large  propeller.  With  a  body  of  aspect  ratio  (length  to  diameter)  2, 
it  is  practically  a  blunt  body  wake.  Figure  20  shows  a  comparison  of  the 
calculation  using 

UN  =  °-3 

Ft  =  2.0 

0  /a  =  0.15 

o 

As  indicated  on  the  figure,  this  initial  turbulent  intensity  level  corresponds 
to  Eto/Edrag  -  1.3/Cp.  Even  with  CD  =  1  for  a  blunt  body  wake,  the  number  is 
still  non-realistic  unless  additional  work  or  energy  input  has  been  required  to 
propell  the  body.  The  possibility  of  additional  drag  on  the  guided  wires  and 
the  non-constant  velocity  nature  of  the  test  can  not  be  discounted.  All  these 
comparisons  tend  to  point  out  that  the  experiment  of  Schooley  and  Stewart, 
though  being  extremely  instructive  in  demonstrating  the  wake  collapse  phenomenon, 
is  quite  far  from  simulating  a  typical  submarine  wake.  For  the  purpose  of 
illustration,  a  calculation,  assuming  the  body  were  streamlined  and  sub-like,  is 
presented  on  the  same  graph.  The  dramatic  difference  in  the  wake  development 
can  explain,  to  a  large  extent,  the  inconsistency  of  the  recent  laboratory 
experiments  with  the  previous  observations  of  Schooley  and  Stewart. 


4.  Summary  and  Concluding  Remarks 


An  improved  phenomenological  model  has  been  developed  for  studying  the 
turbulent  wake  behind  a  submarine  moving  with  a  uniform  speed  in  a  linearly 
stratified  quiet  ocean.  Based  on  a  series  of  parametric  studies,  the  following 
conclusions  are  reached: 

(a)  Initial  anisotropy  has  a  small  effect  on  the  wake  unless  an  unusually 
strong  stratification  is  encountered.  Therefore,  the  initial  turbulent  field 
may  be  assumed  as  isotropic  for  most  practical  situations. 

(b)  Because  of  the  stable  stratification,  the  vertical  turbulence 
component  becomes  vanishingly  small  after  the  wake  reaches  a  maximum  height  and 
the  turbulence  field  becomes  essentially  two-dimensional. 

(c)  The  initial  degree  of  mixing  (gQ)  has  an  appreciable  effect  on  the 
"degree  of  collapse"  as  well  as  the  location  of  collapse  but  it  has  relatively 
small  effect  on  the  other  flow  variables. 

(d)  The  maximum  wake  height  increases  with  increasing  initial  turbulence 
level  UN«  Also,  the  location  of  collapse,  X£  ,  increases  with  increasing  UN  . 

The  choice  of  UN  must  be  consistent  with  the  energy  consideration  and  depends  on 
the  body  configuration. 

(e)  The  effect  of  the.  initial  mean  centerline  velocity  excess  is  quite 

small  on  the  general  flow  field.  A  reasonable  choice  is  U,  =  2  U  . 

do  N 

(f)  Increasing  integral  length  scale  slows  down  the  turbulent  decay  and 
therefore  increases  entrainment.  Since  a  drag  wake  is  expected  to  have  a  larger 
integral  length  scale  than  a  momentumless  wake,  a  larger  b^ay  is  expected  for  a 
drag  wake. 

(g)  Stratification  and  speed  effects  on  the  wake  is  scaled  by  the  Froude 

U0  * 

number  defined  by  Fr  =  —  .  For  increasing  F^,  b^^  increases,  the  physical 
distance  to  collapse  increases  and  the  degree  of  collapse  lessens. 

(h)  A  "pre-collapsed"  wake,  either  caused  by  the  stratification  effect  on 
the  body  boundary  layer  or  by  the  body  geometry,  has  a  strong  effect  on  the  degree 
of  collapse. 

(i)  A  blunt  body  wake  is  expected  to  behave  quite  differently  from  a 
slender  body  sub-like  wake  because  of: 

•  initial  turbulent  intensity 

»  initial  degree  of  mixing 

•  initial  wake  geometry 
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(j)  The  experiment  of  Schooley  and  Stewart  is  atypical  as  compared  to 
a  sub-like  body  wake.  Furthermore,  excessive  drag  or  energy  input  seems  to  be 
deposited  into  the  wake  based  on  the  present  investigation. 

These  conclusions  are  tentative  in  nature,  more  parametric  studies  and 
comparison  with  available  laboratory  experiments  are  required  to  better  assess 
these  conclusions.  The  energy  leakage  through  the  internal  wave  field  must  be 
included  in  the  future.  This  coupling  with  the  internal  wave  field  should  avoid 
the  cvercollapse  phenomenon.  Furthermore,  no  swirl  has  been  included  in  the 
analysis.  In  reality,  a  certain  amount  of  swirl  is  deposited  in  the  wake.  The 
effects  of  swirl  will  have  to  be  investigated  before  the  model  can  be  used  with 
confidence.  A  future  report  will  be  dealing  specifically  with  these  points. 
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Figure  1.  Coordinates  and  Density  Field 
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Figure  7.  Typical  Wake  Area 


Figure  8.  Effect  of  Initial  Degree  of  Mixing  on  the  Wake  Geometry 


Figure  11.  b  v.s.  X/D  for  Various 
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Figure  20.  Comparison  with  the  Experiment  of  Schooley  and  Stewart 


